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ORIGINAL ARTICLE 

A rare WNT1 missense variant overrepresented in ASD leads to 
increased Wnt signal pathway activation 

P-M Martin 1 , X Yang 1 , N Robin 2 , E Lam 3 , JS Rabinowitz 2 , CA Erdman 1 ' 3 , J Quinn 1 ' 3 , LA Weiss 1 ' 3 , SP Hamilton 1,B , P-Y Kwok 3 , RT Moon 2 
and BNR Cheyette 1 

Wnt signaling, which encompasses multiple biochemical pathways that regulate neural development downstream of extracellular 
Wnt glycoprotein ligands, has been suggested to contribute to major psychiatric disorders including autism spectrum disorders 
(ASD). We used next-generation sequencing and Sequenom genotyping technologies to resequence 10 Wnt signaling pathway 
genes in 198 ASD patients and 240 matched controls. Results for single-nucleotide polymorphisms (SNPs) of interest were 
confirmed in a second set of 91 ASD and 1 44 control samples. We found a significantly increased burden of extremely rare missense 
variants predicted to be deleterious by PolyPhen-2, distributed across seven genes in the ASD sample (3.5% in ASD vs 0.8% in 
controls; Fisher's exact test, odds ratio (OR) = 4.37, P = 0.04). We also found a missense variant in WNT1 (S88R) that was 
overrepresented in the ASD sample (8 A/T in 267 ASD (minor allele frequency (MAF) = 1 .69%) vs 1 A/T in 377 controls 
(MAF = 0.13%), OR = 13.0, Fisher's exact test, P = 0.0048; OR = 8.2 and P = 0.053 after correction for population stratification). 
Functional analysis revealed that WNT1-S88R is more active than wild-type WNT1 in assays for the Wnt/p-catenin signaling 
pathway. Our findings of a higher burden in ASD of rare missense variants distributed across 7 of 10 Wnt signaling pathway genes 
tested, and of a functional variant at the WNT1 locus associated with ASD, support that dysfunction of this pathway contributes to 
ASD susceptibility. Given recent findings of common molecular mechanisms in ASD, schizophrenia and affective disorders, these 
loci merit scrutiny in other psychiatric conditions as well. 
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INTRODUCTION 

Autism spectrum disorders (ASD) have a strong genetic etio- 
logy. The last decade of research in psychiatric genetics has led 
to groundbreaking discoveries, including the identification of 
ASD-associated sequence and structural variants in a variety of 
genes implicated in brain development — particularly in synapse 
formation, maintenance and function. However, the mutations 
identified so far explain only a fraction of ASD cases, underscoring 
the need to continue the search for missing genetic links in 
ASD. 

Wnt intercellular signaling refers to a set of key developmental 
pathways that help to organize the nervous system. Among other 
roles, Wnt signaling influences subdivision of the anterior neural 
tube into major brain regions (such as forebrain, midbrain and 
hindbrain), neural precursor proliferation, neural cell fate and 
migration, axon guidance, dendrite development and synapse 
formation. 1 Given this broad spectrum of roles, it is no surprise 
that dysregulation of Wnt signaling causes a similarly broad 
spectrum of deleterious effects on neural development, and 
largely on this basis it has long been considered a candidate 
pathway in psychiatric pathogenesis. Its candidacy in this regard is 
supported by evidence that Wnt signaling pathway components 
are pharmacological targets of lithium, 2,3 hallucinogenic drugs 4 
and antipsychotic medications, 5,6 as well as by evidence derived 
from behavioral tests in animal models 7-11 and from positive 
associations in scattered human genetic studies. 12-17 



Here we used massively parallel sequencing technology to 
simultaneously search for rare variants in 10 Wnt signaling 
pathway genes in a sample of Caucasian patients with ASD and 
healthy controls. We selected these 10 genes out of a much larger 
set of pathway molecules based primarily on empirical support for 
their individual relevance to psychiatry in human genetic {DKK1, 
DKK4, NEUROG1, 07VA/S7), 18-21 animal behavior {DVL1, D/XDC7), 7,10 
biochemical (WNT1, NEUROD1) 22,23 or other types of studies 
{PRICKLE4, DACT2). 24,25 Compared with whole exome sequencing 
strategies that are now underway, this very targeted approach 
benefited from our focus on functionally linked loci for which we 
have substantial expertise, allowing us to immediately follow up 
the most promising variant in in vitro functional assays. 

MATERIALS AND METHODS 

ASD samples (Caucasians) 

Samples were lymphoblastoid cell lines from the Simons Simplex 
Collection. The diagnosis in these samples was based on the ADI-R 
(Autism Diagnostic Interview-Revised). 26 

Unscreened controls (Caucasians) 

Healthy control subjects were recruited by the survey research company 
(Knowledge Networks, Menlo Park, CA, USA) from a nationally representa- 
tive internet-based panel that was selected by random digit dialing. 27 
Subjects completed an online version of the CIDI-SF (Composite 
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International Diagnostic Interview-Short Form) for lifetime history of 
anxiety, mood and substance use disorders. Subjects consented to 
anonymization and deposition of their clinical information and DNA in 
the National Institute of Mental Health (NIMH) repository for use in any 
medical research. 

Pooling samples 

Concentration of individual DNA samples was measured using Pico Green. 
Based on previous reports using an identical pooling strategy, 28,29 we 
expected sequencing error rates to be ~0.5%. To distinguish a single real 
mutant allele in 876 chromosomes (minor allele frequency (MAF) = 0.11%) 
from random sequencing artifacts, we pooled samples into 1 1 sets of 38- 
40 individuals each. Therefore, the frequency of a single mutant allele in 
each pool was 1%, significantly above the detection threshold of 0.5%. 

Target enrichment 

Each pool of genomic DNA from 38-40 individuals was enriched by PCR 
amplification for 10 genes involved in Wnt signaling {WNT1, DVL1, DIXDC1, 
PRICKLE4, DACT2, NEUROG1, DKK1, DKK4, NEUROD1, and CTNNB1), plus 1 
control gene, DACT1, during the validation step (Supplementary Table S1). 
For each gene, the entire locus was amplified in 1 or 2 'long-range' PCR 
reactions. The total amplification length for all 1 1 genes was 71.5 kb. 

Sample preparation and high-throughput sequencing procedure 
Each pool of target-enriched genomic DNA was used to build a prepped 
DNA library following lllumina sample preparation protocols. In brief, each 
pool was ligated with adapters containing a unique index allowing us to 
combine 11 pools in each sequencing lane of the HiSeq flow cell (HiSeq 
2000, lllumina, San Diego, CA, USA). Finished flow cells carrying 438 
individuals (1 lane containing 11 indexed pools of 38-40 individuals) were 
sequenced in the HiSeq 2000 for paired-end analysis. 

Initial data processing (sequence assembly, SNP calling, quality 
assessment) 

Bowtie 30 was used for read alignment against the latest human genome 
reference. SAMtools, 31 Picard (not published; http://picard.sourceforge. 
net/), BEDTools 32 and the Genome Analysis Toolkit (GATK; 33 ) were used for 
post-alignment processing. Multisample realignment around potential 
insertion/deletions (indels) and base quality score recalibration were 
performed before variant calling. For single-nucleotide polymorphism 
(SNP) calling, we used SyZyGy, an algorithm specifically designed to call 
bases in pooled samples. 34 

Analysis of conservation and prediction as deleterious 
For each SNP, analysis by SeattleSeq (http://snp.gs.washington.edu/ 
SeattleSeqAnnotation137/) provided information about base position 
conservation in the genome, whether it was intronic or exonic, and if 
exonic whether it was synonymous, missense or nonsense. Predicted 
deleteriousness for each SNP in exonic sequence was determined by 
PolyPhen-2 35 version v2.1.0r367 or more recent. If either version predicted 
that the SNP was deleterious, it was counted as deleterious. 

Human cDNAs and mutagenesis 

Human WNT1 (clone ID: 30915309) complementary DNA (cDNA) was 
purchased from Thermo Scientific (Rockford, IL, USA) and subcloned into a 
pcDNA 3.1 (-) vector. The WNT1 S88R variant was generated using the 
QuikChange site-directed mutagenesis kit (Stratagene, La Jolla, CA, USA) 
and confirmed by Sanger sequencing. In subsequent functional assays, in 
order to avoid potential confounds arising from nonsystemic biases in 
measurements of plasmid yield or purity, each construct was prepared 
independently >2 times each by different experimenters in two different 
laboratories. Similar functional results were obtained from each indepen- 
dent preparation compared with wild type (WT). 

Luciferase reporter assays 

This was performed essentially as described previously. 36 Briefly, both 
HEK293T and SH-SY5Y BAR-Luc cells were transfected with 1, 5, 20 and 
50 ng pcDNA3.1(-) vector containing either WT WNT1 or WNT1-S88R. 
HEK293T cells were additionally transfected with 20 ng per well pBAR-Luc 
plasmid and 10 ng per well pRL-TK {Renilla luciferase control plasmid). Cells 



were cultured for 24 h, and then luciferase activity measured using the 
Dual-Luciferase Reporter Assay System (Promega, Madison, Wl, USA) and 
an Envision multi-label plate reader (PerkinElmer, Waltham, MA, USA). 

Quantitative reverse transcription-PCR 

Total RNA from cultured HEK293T cells transfected with an empty 
pcDNA3.1(-) vector, or a pcDNA3.1(-) vector containing either WT 
WNT1 or WNT1-S88R, was isolated with the RNeasy Mini kit according to 
the manufacturer's instructions (Qiagen, Hilden, Germany). The cDNA was 
synthesized from 2[ig of total RNA with Superscript III Reverse 
Transcriptase (Invitrogen, Carlsbad, CA, USA) using random hexamer 
primers. Quantitative reverse transcription PCR was performed using the 
PerfeCTa SYBR Green FastMix, ROX (Quanta, Byfleet, UK) with gene-specific 
primer pairs (WNT1 forward (Fwd): S'-GCAGCGACAACATTGACTTC-S', WNT1 
reverse (Rev): S'-GTGGCACTTGCACTCCTG-S'; AXIN2 Fwd: 5'-TTATGCTTTG- 
CACTACGTCCCTCCA-3'; AXIN2 Rev: S'-CGCAACATGGTCAACCCTCAGAC^) 
on an ABI 7900HT real-time PCR system (ABI Advanced Technologies, 
Columbia, MD, USA). Data were analyzed using the comparative CT 
method (User Bulletin No. 2, PerkinElmer Life Sciences, Boston, MA, USA). 
To normalize for loaded cDNA, glyceraldehyde 3-phosphate dehydro- 
genase (GAPDH) was used as an endogenous control {GAPDH Fwd: 
5'<MGGTCATCCATGACAACTTTG-3', GAPDH Rev: 5'-GGCCATCCACAGTCTT 
CTGG-3'). Results shown reflect combined data from two independent 
experiments, each done in triplicate reactions. 

Statistical analyses 

Unless stated otherwise in Figure legends, the data are expressed as 
mean ± s.e.m. Data were analyzed by Student's f-test, Fisher's exact test or 
Pearson's correlation test where appropriate, using GraphPad Prism 
software (GraphPad Software, La Jolla, CA, USA). PLINK software 37 was 
used for Multidimensional Scaling and logistic regression. A P-value of 
<0.05 was considered to be significant. 



RESULTS 

Using a next-generation (massively parallel lllumina-based) 
strategy, we sequenced promoter, 3' untranslated, intronic and 
coding regions of 10 Wnt pathway genes (WNT1, DVL1, DIXDC1, 
PRICKLE4, DACT2, NEUROG1, DKK1, DKK4, NEUROD1 and CTNNB1) in 
genomic DNA obtained from a first set of 198 Caucasian ASD 
patients and 240 healthy Caucasian controls, arranged in pools of 
38-40 samples each. The number of reads was homogeneously 
distributed over the different amplicons and pools (and therefore 
between ASD and controls), reads mapped 88% of the targeted 
regions and coverage averaged 865-fold/individual. To validate 
the next-generation sequencing strategy employed, we compared 
results obtained using this strategy to results obtained using 
Sanger sequencing of the exons in a control locus (DACT1) in a 
subset of samples (40 ASD). Every DACT1 variant {n = 8), including 
several singleton variants (n = 3), found by Sanger sequencing in 
this set was also detected by the next-generation sequencing 
strategy, and the allele frequencies of all SNPs were identical using 
each methodology (Pearson's correlation test; r = 0.99; P< 0.0001; 
Supplementary Table S1). 

In the 10 genes under study, we identified 652 SNPs overall in 
this set of 438 combined ASD and control samples. We restricted 
most of our attention to nonsynonymous sequence variants found 
only in our own sample set (not referenced in dbSNP, 1000 
Genomes, the Exome Variant Server (EVS) or the Exome Chip 
Design Consortium). To determine whether any of the putative 
SNPs found by this method were artifacts produced by the next- 
generation sequencing technology, we confirmed each of them 
by individual genotyping using either IPlex Sequenom or Sanger 
sequencing. Using this approach we confirmed 12 novel singleton 
missense variants distributed across 8 of the 10 loci included in 
our study (Supplementarty Table S2). Moreover, genotyping in a 
second set of 91 ASD and 144 control samples failed to identify 
any additional occurrences of each of these rare deleterious 
missense SNPs. Although each of these SNPs were found only 
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Figure 1. Burden of novel unique variants in the Wnt pathway in 
autism spectrum disorders (ASD). Percentage of individuals carrying 
a novel unique variant that was nonsynonymous and predicted 
to be deleterious by the PolyPhen-2 algorithm, compared 
with nonsynonymous, synonymous, untranslated region (UTR) or 
intronic single-nucleotide polymorphisms (SNPs). *P<0.05; NS, not 
significant. 



once among the probands in our samples and had not been 
reported in any public database, genotyping of parents demon- 
strated that none of them had occurred de novo; all were 
transmitted from either the mother or father. Moreover, sequen- 
cing of family trios showed that each was also present in an 
unaffected sibling. This verifies that each of these variants is not 
an idiosyncratic artifact arising during lymphoblastoid cell line 
derivation, and suggests that they are extremely rare alleles 
(single-nucleotide variants) in the Caucasian population repre- 
sented by our samples. 

Although we found a greater number of extremely rare 
missense variants in our ASD sample compared with our control 
sample (4% in ASD vs 1.7% in controls; Figure 1), this difference in 
the 'rare missense variant burden' did not reach statistical 
significance (Fisher's exact test, P = 0.1 1). Nonetheless, when 
considering only those rare missense variants that were also 
predicted to be deleterious by PolyPhen-2, the difference in ASD 
versus control groups did become significant (3.5% in ASD vs 0.8% 
in controls; Fisher's exact test, OR = 4.37, P = 0.04; Figure 1). We 
also examined the distribution of less informative variants (that is, 
synonymous, untranslated regions and intronic variants) in each 
group, and found a very similar distribution for all of them 
(synonymous: 3.5% in ASD vs 2.9% in controls; untranslated 
regions: 1 1.1% in ASD vs 9.6% in controls; intronic: 24% in ASD vs 
21% in controls; Figure 1). The similarity in the incidence of less 
informative variants argues against the possibility that the 
difference in the incidence of rare deleterious variants can be 
attributed to differences in recent ancestry between the ASD and 
control sample groups. 

In addition to these novel SNPs, we found one previously 
referenced rare (<1%) missense SNP that was overrepresented in 
our ASD group versus our control group. This SNP (rs61 758378) 
occurred in WNT1, the gene encoding the first-discovered 
mammalian Wnt ligand. To eliminate the possibility of a 
sequencing artifact we confirmed this result by genotyping all 
samples using IPlex Sequenom. We also genotyped an additional 
91 ASD and 144 control samples to replicate and potentially 
extend these findings. In this way, we not only confirmed MAF 
results for rs61 758378 obtained using the lllumina platform in the 
original sample set, we further found that the rs61 758378 minor 
allele remained significantly overrepresented in the combined 
ASD sample compared with the combined control sample (8 A/T 
in 267 ASD (MAF = 1 .69%) vs 1 A/T in 377 controls (MAF = 0.1 3%); 



Fisher's exact test, allelic P = 0.0048; Table 1). Furthermore, when 
we compared the rs61 758378 MAF in our combined ASD sample 
with the reported MAF in the EVS for a much larger group of 
individuals, the association also remained significant (8 A/T in 267 
in our combined ASD sample (MAF = 1.69%) vs 45 A/T in 4300 
European/American individuals in the EVS (MAF = 0.5%), Fisher's 
exact test, allelic P = 0.011). Because different allele frequencies 
are reported in European populations (TSI, Tuscans in Italy, 
MAF = 3.0%; CEU, Northern Europeans in Utah, MAF = 0.5%), we 
wanted to evaluate and address possible effects of subtle 
stratification in ancestry between cases and controls. We therefore 
performed Multidimensional Scaling using genome-wide SNP 
genotype data for a subset of the total samples used in this 
analysis (235 SSCs and 329 controls) using PLINK. Multidimen- 
sional Scaling using reference CEU and TSI data showed that the 
ratio of individuals clustering with CEU and TSI was not 
significantly different between ASD and controls. Nevertheless, 
we noted that one of the first four dimensions was significantly 
different between ASD and controls, indicating possible stratifica- 
tion. We sought to account for this using a logistic regression 
analysis including the first four dimensions as covariates. In this 
analysis, the odds ratio (OR) remained high (OR = 8.1 7) and the 
P-value remained close to significance (P = 0.053). We note that 
our analyses have low power for SNPs with such a low MAF and in 
a modest number of samples. 

The rs61 758378 minor allele changes a conserved serine at 
position 88 into arginine. Although this change is not predicted to 
be deleterious by PolyPhen-2 (HumDiv score = 0.1 74), serine 88 is 
nonetheless a very highly conserved residue across species 
(PhasCons score = 0.998), and the change to arginine at this 
residue is predicted to be deleterious by Grantham score, which 
categorizes codon replacements into classes of increasing 
chemical dissimilarity (Grantham score = 110, moderately radical 
change). 38 Given the potential association between the WNT1 
rs61 758378 minor allele and ASD in our sample set, we therefore 
tested for functional effects of this variant on Wnt signaling 
pathway activation. To accomplish this, we engineered this variant 
into a human WNT1 cDNA (WNT1-S88R) using site-directed 
mutagenesis and then tested it alongside WT Wnt1 in several 
assays for Wnt signaling activity. 

First, we performed a plasmid-based Wnt reporter assay in 
human embryonic kidney (HEK293T) tissue culture cells. This is a 
well-established method to quantify transcriptional activation 
downstream of the Wnt/(3-catenin signaling pathway, based in this 
case on transient transfection of the pBAR-Luc reporter plasmid 
that expresses luciferase downstream of a (3-catenin-responsive 
promoter. 36 Cells were co-transfected with pBAR-Luc and a dose 
range (1 -50 ng) of plasmids encoding WNT1 -S88R or WT WNT1 to 
test for the ability of these proteins to activate Wnt/(3-catenin 
pathway-dependent transcription. At all doses tested, WNT1-S88R 
showed increased activity in this assay compared with WT WNT1 
(Figure 2a). Moreover, this difference was statistically significant at 
the 20 ng dose (WT WNT1 , 1 .00 ± 0.09 vs WNT1 -S88R, 1.37 ± 0.087, 
f-test, P=0.02) and at the 50 ng dose (WT WNT1, 1.00 ±0.18 vs 
WNT1-S88R, 1.86 ±0.22, f-test, P = 0.016; Figure 2a). Similar results 
were obtained when experiments were replicated in a different 
laboratory by a different investigator using independently pre- 
pared DNA constructs and reagents (pSuperTop; data not shown). 

The aforementioned assays rely on transient transfection of 
reporter plasmid, and on this basis they are conceptually more 
variable than assays based on read-out of an endogenous Wnt 
pathway target gene or of a single-copy Wnt pathway reporter 
integrated into the genome. Accordingly, we confirmed these 
results in the same (HEK293T) cell line using quantitative reverse 
transcription-PCR to directly measure transcriptional activation of 
the endogenous Wnt/(3-catenin pathway target gene AXIN2. 39 
Mirroring results obtained with the Wnt reporter assay, com- 
pared with WT WNT1, WNT1-S88R showed significantly increased 
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Table 1. Association of WNT1 rs61 758378 with ASD 
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Figure 2. Functional analysis of S88R. (a) Plasmid-based Wnt/p- 
catenin pathway reporter assay in HEK293T cells. HEK293T cells were 
transfected with the pBAR-Luc reporter plasmid together with 
plasmids encoding either WNT1-S88R or wild-type (WT) WNT1 at 
doses ranging from 1 to 50 ng. (b) Quantitative reverse transcription- 
PCR (RT-PCR) of the Wnt/p-catenin pathway endogenous target 
gene AXIN2 in HEK293T cells transfected with 5 ng of each WNT1 
plasmid, expressed as a ratio to WNT1 expression also determined 
by RT-PCR. (c) Luciferase reporter activation in SH-SY5Y BAR-Luc 
cells transfected with 50 ng of each WNT1 plasmid. *P<0.05 and 
***P< 0.001. 



transcriptional activation of AXIN2 (WT Wnt1, 1 .00 ± 0.025 vs Wnt1- 
S88R, 2.1610.43, f-test, P = 0.02; Figure 2b). We confirmed these 
results via a third in vitro assay that relies on SH-SY5Y BAR-Luc 



cells. This is an independently derived immortalized cell line with a 
single copy of a Wnt/(3-catenin pathway luciferase reporter stably 
integrated into its genome. 36 Consistent with data obtained by 
the plasmid-based Wnt reporter assays and by quantitative 
reverse transcription-PCR in HEK293T cells, SH-SY5Y BAR-Luc 
cells showed significantly increased activation of their reporter 
after transfection of WNT1-S88R (WT Wntl, 1.00 ±0.04 vs Wnt1- 
S88R, 1.28 ±0.03, f-test, P = 0.0003; Figure 2c). 

Taken together, our data from cell-based transcription target 
assays suggest that the rs61 758378 minor allele associated with 
ASD in our sample (WNT1-S88R) encodes a modestly more active 
form of the WNT1 ligand compared with the most prevalent WNT1 
allele found in Caucasians. 



DISCUSSION 

Our finding in ASD of a higher burden of rare missense variants 
distributed across 7 of 10 Wnt signaling pathway genes tested, 
and of an ASD-associated functional variant at the WNT1 locus, 
supports that dysfunction of this pathway contributes to ASD 
susceptibility. 

In functional assays for the Wnt/(3-catenin pathway, the 
biochemical cascade physiologically activated by the WNT1 gene 
product, 40,41 WNT1-S88R consistently showed increased activity. In 
animal models, hyperactivation of Wnt/(3-catenin signaling can 
lead to dominant neurodevelopmental defects including in 
brain regionalization 42,43 neural precursor proliferation, cell fate 
specification, neuron migration, 44-47 and in dendrite and synapse 
formation or function. 48,49 A mutation in a gene, such as WNT1, that 
modestly increases activation of this downstream signaling pathway 
would therefore be expected to result in neurodevelopmental and 
behavioral phenotypes. Moreover, as ASD and other major 
psychiatric disorders are generally not associated with gross 
neuroanatomical defects, changes in Wnt signaling activity 
contributing to susceptibility for these disorders should be modest. 
Anatomical changes resulting from a signaling variant such as WNT1- 
S88R might only be evident at a cellular level (for example, in neural 
subtype numbers and position, neural morphology, connectivity and 
so on) or at a subcellular level (for example, in dendritic spine 
morphology, synapse numbers, synapse strength, synapse plasticity 
and so on), yet could still be sufficient to cause behavioral symptoms. 
This of course conforms closely to the central features of major 
psychiatric disorders, which are distinguished from neurological 
disorders in part by the absence of focal or grossly evident structural 
brain defects. Moreover, it is not our intent to posit that this or the 
other rare coding variants we have identified here in WNT1 and in 
Wnt signaling pathway component genes are monogenic causes of 
psychiatric illness. 50 Genetic susceptibility to these disorders is highly 
heterogeneous and thought to be frequently multigenic; that is, for 
most patients genetic loading is thought to be distributed across 
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many loci that affect a common pathogenic pathway. 51,52 In this 
susceptibility model, an allelic variant at any one locus may only 
appreciably increase disease risk in the context of variants at other 
loci that affect the same pathway in the same individual. The rare 
nonsynonymous sequence variants found in both affected and 
unaffected individuals within the same family in this study may 
contribute to ASD in the context of other genetic and nongenetic 
factors that increase susceptibility. 

Several lines of evidence support that Wnt ligands (of which 
there are 19 family members in humans) participate in ASD 
susceptibility. A recent whole exome sequencing study identified 
copy number variants affecting both WNT3 and WNT9B in an ASD 
patient, 53 whereas another study showed that WNT3 expression is 
increased in the prefrontal cortex of young autistic patients 
compared with controls. 54 Rare variants in WNT2 have been 
reported to segregate with ASD. 16 Although other groups have 
not always replicated this linkage, 55,56 two more recent studies 
have found SNPs in WNT2 associated with either ASD or speech 
delay. 15,17 In a Valproate-exposure rat model of ASD, expression of 
both WNT1 and WNT2 was upregulated in the developing 
forebrain and this was associated with increased Wnt/(3-catenin 
pathway activity in that brain region. 57 Other recent sequencing 
studies have identified several ASD-susceptibility genes that 
encode functionally linked components of a protein network 
involved in the regulation of (3-catenin-dependent transcription, 
as occurs in Wnt/(3-catenin signaling. 21,58,59 The findings of the 
present study add to this evidence by suggesting that minor 
alleles and extremely rare variants at the WNT1 locus and other 
Wnt pathway loci may contribute to risk for this psychiatric 
disorder, a possibility that merits further attention in ongoing 
human genetic studies focused on ASD. Given evidence that 
genetic etiologies are shared across several psychiatric 
disorders, 60-62 as well as that Wnt signaling is involved more 
broadly in psychiatric pathogenic and therapeutic mechanisms, 63 
these loci should also be carefully scrutinized in ongoing human 
genetic studies of other major psychiatric disorders, such as 
schizophrenia and major affective disorders. 



CONFLICT OF INTEREST 

The authors declare no conflict of interest. 



ACKNOWLEDGEMENTS 

P-MM has received support from a Pilot Project award sponsored by Kristine Yaffe, 
the Roy and Marie Scola Endowed Chair in Psychiatry at UCSF; XYY was funded by 
NIH T32 MH089920; this research has also been supported by grants to BNRC and 
SPH together and to BNRC alone from the UCSF Research Allocation Program (RAP) 
including the Program for Breakthrough Biomedical Research (PBBR), by the 
Department of Psychiatry, and by the Center for Neurobiology and Psychiatry. LAW 
has support from the International Mental Health Research Organization. We 
acknowledge SFARI for ASD genomic DNA samples and for phenotypic data made 
available on SFARI Base. We also thank all the families who donated samples and the 
principal investigators involved in their collection (A Beaudet, R Bernier, 
J Constantino, E Cook, E Fombonne, D Geschwind, D Grice, A Klin, D Ledbetter, 
C Lord, C Martin, D Martin, R Maxim, J Miles, O Ousley, B Peterson, J Piggot, C Saulnier, 
M State, W Stone, J Sutcliffe, C Walsh and E Wijsman). Approved researchers can 
obtain the SSC population data set by applying at https://base.sfari.org. Collection of 
the control data set was supported by grants from the NIMH and the National 
Alliance for Research on Schizophrenia and Depression. 



REFERENCES 

1 Mulligan KA, Cheyette BN. Wnt signaling in vertebrate neural development and 
function. J Neuroimmune Pharmacol 2012; 7: 774-787. 

2 Wexler EM, Geschwind DH, Palmer TD. Lithium regulates adult hippocampal 
progenitor development through canonical Wnt pathway activation. Mol Psy- 
chiatry 2008; 13: 285-292. 

3 O'Brien WT, Klein PS. Validating GSK3 as an in vivo target of lithium action. 
Biochem Soc Trans 2009; 37: 1 133-1 138. 



4 Svenningsson P, Tzavara ET, Carruthers R, Rachleff I, Wattler S, Nehls M et al. 
Diverse psychotomimetics act through a common signaling pathway. Science 
2003; 302: 1412-1415. 

5 Freyberg Z, Ferrando SJ, Javitch JA. Roles of the Akt/GSK-3 and Wnt signaling 
pathways in schizophrenia and antipsychotic drug action. Am J Psychiatry 2010; 
167: 388-396. 

6 Sutton LP, Honardoust D, Mouyal J, Rajakumar N, Rushlow W. Activation of the 
canonical Wnt pathway by the antipsychotics haloperidol and clozapine involves 
dishevelled-3. J Neurochem 2007; 102: 153-169. 

7 Long JM, LaPorte P, Paylor R, Wynshaw-Boris A. Expanded characterization of the 
social interaction abnormalities in mice lacking Dvll. Genes Brain Behav 2004; 3: 
51-62. 

8 Gould TD, O'Donnell KC, Picchini AM, Dow ER, Chen G, Manji HK. Generation and 
behavioral characterization of beta-catenin forebrain-specific conditional knock- 
out mice. Behav Brain Res 2008; 189: 117-125. 

9 Mines MA, Yuskaitis CJ, King MK, Beurel E, Jope RS. GSK3 influences social pre- 
ference and anxiety-related behaviors during social interaction in a mouse model 
of fragile X syndrome and autism. PLoS One 2010; 5: e9706. 

10 Kivimae S, Martin PM, Kapfhamer D, Ruan Y, Heberlein U, Rubenstein JL et al. 
Abnormal behavior in mice mutant for the Disci binding partner, Dixdd. Transl 
Psychiatry 2011; 1: e43. 

1 1 Sowers LP, Loo L, Wu Y, Campbell E, Ulrich JD, Wu S et al. Disruption of the non- 
canonical Wnt gene PRICKLE2 leads to autism-like behaviors with evidence for 
hippocampal synaptic dysfunction. Mol Psychiatry advance online publication, 
28 May 2013; doi:1 0.1 038/mp.201 3.71 (e-pub ahead of print). 

12 Cui DH, Jiang KD, Jiang SD, Xu YF, Yao H. The tumor suppressor adenomatous 
polyposis coli gene is associated with susceptibility to schizophrenia. Mol Psy- 
chiatry 2005; 10: 669-677. 

13 Mao Y, Ge X, Frank CL, Madison JM, Koehler AN, Doud MK et al. Disrupted in 
schizophrenia 1 regulates neuronal progenitor proliferation via modulation of 
GSK3beta/beta-catenin signaling. Cell 2009; 136: 1017-1031. 

14 Katsu T, Ujike H, Nakano T, Tanaka Y, Nomura A, Nakata K et al. The human 
frizzled-3 (FZD3) gene on chromosome 8p21, a receptor gene for Wnt ligands, is 
associated with the susceptibility to schizophrenia. Neurosci Lett 2003; 353: 53-56. 

15 Lin PI, Chien YL, Wu YY, Chen CH, Gau SS, Huang YS et al. The WNT2 gene 
polymorphism associated with speech delay inherent to autism. Res Dev Disabil 
2012 Res Dev Disabil 33: 1533-1540. 

16 Wassink TH, Piven J, Vieland VJ, Huang J, Swiderski RE, Pietila J et al. Evidence 
supporting WNT2 as an autism susceptibility gene. Am J Med Genet 2001; 105: 
406-41 3. 

17 Marui T, Funatogawa I, Koishi S, Yamamoto K, Matsumoto H, Hashimoto O et al. 
Association between autism and variants in the wingless-type MMTV integration 
site family member 2 ( WNT2) gene. Int J Neuropsychopharmacol 2010; 13: 
443-449. 

18 Aleksic B, Kushima I, Ito Y, Nakamura Y, Ujike H, Suzuki M et al. Genetic association 
study of KREMEN1 and DKK1 and schizophrenia in a Japanese population. Schi- 
zophr Res 201 0; 1 1 8: 1 1 3-1 1 7. 

19 Fanous AH, Chen X, Wang X, Amdur RL, O'Neill FA, Walsh D et al. Association 
between the 5q31.1 gene neurogeninl and schizophrenia. Am J Med Genet B 
Neuropsychiatr Genet 2007; 144B: 207-214. 

20 Proitsi P, Li T, Hamilton G, Di Forti M, Collier D, Killick R et al. Positional pathway 
screen of wnt signaling genes in schizophrenia: association with DKK4. Biol Psy- 
chiatry 2008; 63: 13-16. 

21 O'Roak BJ, Vives L, Girirajan S, Karakoc E, Krumm N, Coe BP et al. Sporadic autism 
exomes reveal a highly interconnected protein network of de novo mutations. 
Nature 2012; 485: 246-250. 

22 Miyaoka T, Seno H, Ishino H. Increased expression of Wnt-1 in schizophrenic 
brains. Schizophr Res 1999; 38: 1-6. 

23 Aston C, Jiang L, Sokolov BP. Microarray analysis of postmortem temporal cortex 
from patients with schizophrenia. J Neurosci Res 2004; 77: 858-866. 

24 Fisher DA, Kivimae S, Hoshino J, Suriben R, Martin PM, Baxter N et al. Three Dact 
gene family members are expressed during embryonic development and in the 
adult brains of mice. Dev Dyn 2006; 235: 2620-2630. 

25 Veeman MT, Slusarski DC, Kaykas A, Louie SH, Moon RT. Zebrafish prickle, a 
modulator of noncanonical Wnt/Fz signaling, regulates gastrulation movements. 
CurrBiol 2003; 13: 680-685. 

26 Lord C, Rutter M, Le Couteur A. Autism Diagnostic Interview-Revised: a revised 
version of a diagnostic interview for caregivers of individuals with possible per- 
vasive developmental disorders. J Autism Dev Disord 1994; 24: 659-685. 

27 Sanders AR, Levinson DF, Duan J, Dennis JM, Li R, Kendler KS et al. The Internet- 
based MGS2 control sample: self report of mental illness. Am J Psychiatry 2010; 
167: 854-865. 

28 Calvo SE, Tucker EJ, Compton AG, Kirby DM, Crawford G, Burtt NP et al. High- 
throughput, pooled sequencing identifies mutations in NUBPL and FOXRED1 in 
human complex I deficiency. Nat Genet 2010; 42: 851-858. 



© 2013 Macmillan Publishers Limited 



Translational Psychiatry (2013), 1 -6 



Functional WNT1 variant in autism 
P-M Martin et a I 



29 Nejentsev S, Walker N, Riches D, Egholm M, Todd JA. Rare variants of IFIH1, a gene 
implicated in antiviral responses, protect against type 1 diabetes. Science 2009; 
324: 387-389. 

30 Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient align- 
ment of short DNA sequences to the human genome. Genome Biol 2009; 10: R25. 

31 Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N et al. The Sequence 
Alignment/Map format and SAMtools. Bioinformatics 2009; 25: 2078-2079. 

32 Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic 
features. Bioinformatics 2010; 26: 841-842. 

33 DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C et al. A frame- 
work for variation discovery and genotyping using next-generation DNA 
sequencing data. Nat Genet 2011; 43: 491-498. 

34 Rivas MA, Beaudoin M, Gardet A, Stevens C, Sharma Y, Zhang CK et al. Deep 
resequencing of GWAS loci identifies independent rare variants associated with 
inflammatory bowel disease. Nat Genet 2011; 43: 1066-1073. 

35 Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P et al. 
A method and server for predicting damaging missense mutations. Nat Methods 
2010; 7: 248-249. 

36 Biechele TL, Moon RT. Assaying beta-catenin/TCF transcription with beta-catenin/ 
TCF transcription-based reporter constructs. Methods Mol Biol 2008; 468: 99-1 10. 

37 Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D et al. PLINK: a 
tool set for whole-genome association and population-based linkage analyses. 
Am J Hum Genet 2007; 81: 559-575. 

38 Grantham R. Amino acid difference formula to help explain protein evolution. 
Science 1974; 185: 862-864. 

39 Jho EH, Zhang T, Domon C, Joo CK, Freund JN, Costantini F. Wnt/beta-catenin/Tcf 
signaling induces the transcription of Axin2, a negative regulator of the signaling 
pathway. Mol Cell Biol 2002; 22: 1 172-1 183. 

40 Moon RT, Brown JD, Torres M. WNTs modulate cell fate and behavior during 
vertebrate development. Trends Genet 1997; 13: 157-162. 

41 Willert K, Nusse R. Beta-catenin: a key mediator of Wnt signaling. Curr Opin Genet 
Dev 1998; 8: 95-102. 

42 Fekany-Lee K, Gonzalez E, Miller-Bertoglio V, Solnica-Krezel L. The homeobox gene 
bozozok promotes anterior neuroectoderm formation in zebrafish through negative 
regulation of BMP2/4 and Wnt pathways. Development 2000; 127: 2333-2345. 

43 Takai A, Inomata H, Arakawa A, Yakura R, Matsuo-Takasaki M, Sasai Y. Anterior neural 
development requires Dell, a matrix-associated protein that attenuates canonical 
Wnt signaling via the Ror2 pathway. Development 2010; 137: 3293-3302. 

44 Tang K, Yang J, Gao X, Wang C, Liu L, Kitani H et al. Wnt-1 promotes neuronal 
differentiation and inhibits gliogenesis in P19 cells. Biochem Biophys Res Commun 
2002; 293: 167-173. 

45 Ligon KL, Echelard Y, Assimacopoulos S, Danielian PS, Kaing S, Grove EA et al. Loss 
of Emx2 funtion leads to ectopic expression of Wnt1 in the developing tele- 
ncephalon and cortical dysplasia. Development 2003; 130: 2275-2287. 

46 Bhat KM. Wingless activity in the precursor cells specifies neuronal migratory 
behavior in the Drosophila nerve cord. Dev Biol 2007; 311: 613-622. 

47 Joksimovic M, Patel M, Taketo MM, Johnson R, Awatramani R. Ectopic Wnt/beta- 
catenin signaling induces neurogenesis in the spinal cord and hindbrain floor 
plate. PLoS One 2012; 7: e30266. 



48 Yu X, Malenka RC. Beta-catenin is critical for dendritic morphogenesis. Nat 
Neurosci 2003; 6: 1 1 69-1 1 77. 

49 Ataman B, Ashley J, Gorczyca M, Ramachandran P, Fouquet W, Sigrist SJ et al. 
Rapid activity-dependent modifications in synaptic structure and function require 
bidirectional Wnt signaling. Neuron 2008; 57: 705-718. 

50 De Ferrari GV, Moon RT. The ups and downs of Wnt signaling in prevalent neu- 
rological disorders. Oncogene 2006; 25: 7545-7553. 

51 Buxbaum JD, Daly MJ, Devlin B, Lehner T, Roeder K, State MW. The autism 
sequencing consortium: large-scale, high-throughput sequencing in autism 
spectrum disorders. Neuron 2012; 76: 1052-1056. 

52 Devlin B, Scherer SW. Genetic architecture in autism spectrum disorder. Curr Opin 
Genet Dev 2012; 22: 229-237. 

53 Levy D, Ronemus M, Yamrom B, Lee YH, Leotta A, Kendall J et al. Rare de novo and 
transmitted copy-number variation in autistic spectrum disorders. Neuron 2011; 
70: 886-897. 

54 Chow ML, Pramparo T, Winn ME, Barnes CC, Li HR, Weiss L et al. Age-dependent 
brain gene expression and copy number anomalies in autism suggest distinct 
pathological processes at young versus mature ages. PLoS Genet 2012; 8: 
e1 002592. 

55 McCoy PA, Shao Y, Wolpert CM, Donnelly SL, Ashley-Koch A, Abel HL et al. No 
association between the WNT2 gene and autistic disorder. Am J Med Genet 2002; 
114: 106-109. 

56 Li J, Nguyen L, Gleason C, Lotspeich L, Spiker D, Risch N et al. Lack of evidence for 
an association between WNT2 and RELN polymorphisms and autism. Am J Med 
Genet B Neuropsychiatr Genet 2004; 126B: 51-57. 

57 Wang Z, Xu L, Zhu X, Cui W, Sun Y, Nishijo H et al. Demethylation of specific Wnt/ 
beta-catenin pathway genes and its upregulation in rat brain induced by prenatal 
valproate exposure. Anat Rec (Hoboken) 2010; 293: 1947-1953. 

58 O'Roak BJ, Vives L, Fu W, Egertson JD, Stanaway IB, Phelps IG et al. Multiplex 
targeted sequencing identifies recurrently mutated genes in autism spectrum 
disorders. Science 2012; 338: 1619-1622. 

59 Neale BM, Kou Y, Liu L, Ma'ayan A, Samocha KE, Sabo A et al. Patterns and rates 
of exonic de novo mutations in autism spectrum disorders. Nature 2012; 485: 
242-245. 

60 State MW, Levitt P. The conundrums of understanding genetic risks for autism 
spectrum disorders. Nat Neurosci 2011; 14: 1499-1506. 

61 Talkowski ME, Rosenfeld JA, Blumenthal I, Pillalamarri V, Chiang C, 
Heilbut A et al. Sequencing chromosomal abnormalities reveals neurodeve- 
lopmental loci that confer risk across diagnostic boundaries. Cell 2012; 149: 
525-537. 

62 Smoller JW. Identification of risk loci with shared effects on five major psychiatric 
disorders: a genome-wide analysis. Lancet 2013; 381: 1371-1379. 

63 Okerlund ND, Cheyette BN. Synaptic Wnt signaling-a contributor to major psy- 
chiatric disorders? J Neurodev Disord 201 1; 3: 162-174. 



This work is licensed under a Creative Commons Attribution- 
NonCommercial-NoDerivs 3.0 Unported License. To view a copy of 



this license, visit http://creativecommons.Org/licenses/by-nc-nd/3.0/ 



Supplementary Information accompanies the paper on the Translational Psychiatry website (http://www.nature.com/tp) 



Translational Psychiatry (2013), 1 -6 



© 2013 Macmillan Publishers Limited 



